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ABSTRACT 


The mathematical model of infectious diseases is a tool that has been used 
to study the mechanisms by which disease spread, to predict the future and 
course of an outbreak and to evaluate strategies to control an epidemic. We 
envisaged a community of n individuals comprising at time t, x, susceptible, 
y infectious in circulation and z individuals who were isolated, dead, or 
recovered and immune. We further postulated infection and removal rates 
B and y, so that there wave Bxydt new infections and yydt removals in time 
tithe simplest way to do this is to introduce a birth parameter u, so at to 
give dt new susceptible in time dt. If the population is to remain stable the 
arrival of new susceptible must be balanced by an appropriately defined 
birth rate. The present paper represents the model in special way, in which 
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the infection occurs in human's body then the resistance of body gradually 


decays same as motion decays in damped oscillation. On solving the 
equation of model, we get solution that gives the idea about the seasonal 


variation in infection. 
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INTRODUCTION: 

The mathematical model of infectious diseases is a tool 
that has been used to study the mechanisms by which 
disease spread, to predict the future course of an outbreak 
and to evaluate strategies to control an epidemic(Eisinger 
and Thulke, 2008;Parshani, et. al., 2010;Valdez, et. al., 
2020;).We envisaged a community of n individuals 
comprising at time t, x, susceptible, y infectious in 
circulation and z individuals who were isolated, dead, or 
recovered and immune (Bailey, 1964; 1975). We further 
postulated infection and removal rates Pandy, so that 
there wave fxydt new infections and yydt removals in 
time dt. An additional assumption to be made here is that 
the stock of susceptible is continually replenished 
(Barlett, 1956 and Brown, 1971). 


The simplest way to do this is to introduce a birth 
parametery, so at to give dt new susceptible in time dt. If 
the population is to remain stable the arrival of new 
susceptible must be balanced by an appropriately defined 
birthrate(Carrier, 1996; Cosma 2018; Sridhar and 
Majumder, 2020). The simplest model that can be 
constructed avoids explicit reference to a death rate by 
concentrating on the groups of susceptible and infective, 
the former at any rate being supposed not subject to death 
(Bodmer and Cavalli, 1976). This is equivalent to 
assuming that on average the deaths of removed 
individuals are just balanced by births of new susceptible 
(Maki, and Thompson, 1973). 
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MODEL DEVELOPMENT AND APPLICATION: 


dx 

ar BXY +B 

dy _ a 

aT BX YY (1.1) 
dz 

Fee Al 


Equilibrium values Xo and Yo are given by equating the 
differential coefficients to zero. 


Xo 4 
Thus yo= 5 (1.2) 


The equations for small departures from these equilibrium 
values are obtained by writing 


X= Xo (1+ u) (1.3) 
y= yo (1+v) 

and substituting (1.3) in (1.1) 

goose (u+v-+tuv) (1.4) 


dt 


dv 
TET ud ty) 


Where 
¥ 7 
opt ae 


It we now neglect the product uv and eliminate u from the 
two equations (1.4) we obtain the second order 
differential equations in v. 
d?v 1dv 1 

+ — 


+—v=0 (1.5) 


dt? odt oT 
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1 1 
D¢2v+—Dv+—v=0 
oO oT 


or, (D2? ++ =D + —v=0 
o OT 
jie (22* 
_ GF Va oY 
So, D= - 
oer a ee = 
20 2\ oT o 


Solution of the equation (1.5) is 


= ert {C,cos ( | — 4) t+ Csin( |*— +y¢ 
v = ezat {Cicos ( = =e) 2sin( = salt} 


with help of boundary condition for t = 40 


Zz 4 1 
we have v = voe2e cos|_|—— — Jt 
oT oO 


-1 
V1 = Voe2e cosét (1.6) 


Where,é = ( — — 3) 


o2 


for a suitably chosen origin of time, we then obtain the 
solution for u given by 


ui = Vo () ea cos (ét+w) 


1 
T 


Where, cos p = -2 (2)? Os yn (1.7) 


oO 


The solutions wu; and vi, which are linearized or first order 
components only, clearly involve damped harmonic trains 
of waves with period 2 2 /€. In this application to measles, 
T equal to two weeks the approximate incubation period 
and, from the data available to him he estimated o in 
London to be roughly 68.2 weeks. Equation (1.6) give the 
period 2 2 /& = 75.7 weeks, with a peak to peak damping 


Tw 1 
factor of e o €2 =0.58. 


For somewhat larger oscillations we ought to take in to 
account the non-linear character of the equations (1.4). 
For this purpose it is convenient to write: 


u= Uytaq1U12+a12U,U1+a22U712+ aveudstvaiwant 
u= Ui+a41U12+a12U,U1+a22U12+ eusbaeie vty (1.8) 


Where the coefficient aj and bj, can be found in terms of = 
by straight forward substitution in (1.4). It = is small, so 
that w - 2/2 we have approximately: 


u=u'1 (1 + =0;) 
1 =e 
v=vit 3 V70e20 cos2ét (1.9) 
T ies 
ui= -vo(=)'/2e%0 sinét 


It will be noticed that the damping coefficient has 
relatively little influence on the period = which for small < 


1 
is roughly 27(oT)2 = 2 1./ Bu. 


And so largely depends on the birth rate for new 
susceptible and the infection rate. 


numerical solution of (1.1) through special care in needed 
to ensure sufficient accuracy and to avoid mistaken 
conclusion about the effect of damping. 


Another way of representing the oscillatory behavior of 
the process is to plot the path traced by the point (x,y). 
This also allows a convenient method of comparing the 
deterministic solution with path followed by actual 
realizations of the stochastic analogue. It can be shown 
that the deterministic curve found in this way has the form 
of a spiral converging on the equilibrium point (xo,yo), and 
that this occurs is spite of the non-linear form of (1.1). 


An elegant argument due to G.E.H. Reuter (Bortlett, 1956) 
is as follows. 


Consider the function 





F (uv) = {(1+u)-log(1+u)+(=){(1+v)-log(1+v)} (1.10) 
Differentiating with respect to t and using (4) 

af\ _ u? 
(4) = o (1+u) <0 (1.11) 


Thus f continually decreases along any path for which t 
increases. Since i=+ a if follows that f tends to a finite 


limite for fox 1 += as t &. Now the curve f = c are closed 


surround the point (xo,yo) and shrink down as c 1 + = 


By considering the second different coefficient d#f/dtz we 
can show that fo=1+, so that the point (x, y) must actually 
tends to (xo,Yo). 

An important consequence of the above discussion is that, 
while the additional of a constant influx of fresh 
susceptible is sufficient to account for epidemic waves a 
period of about the order of magnitude the damping down 
to a steady endemic state entailed by the calculations is at 
variance with observed epidemiological facts. 


ESTIMATION OF SEASONAL VARIATIONS IN INFECTION 
RATE: 

The simplest modification is to try replacing the infection 
rate 8 by B' = B+ Bi cos wt. 


Where, 27 /w=52 weeks. 


We suppose the relative amplitude of these forced 
oscillations. p = a to be small. 


Substituting the new value of f in (1.1) and using (1.3) 
gives equations corresponding to(1.4). These can now be 
deal with as before by eliminationuto yield to 
modified form of (1.5). 

d?v_ 1dv 1 


TW. 
+ Vv=- sinwt 
dt? odt ot T 





(1.12) 


The particular integral of (1.12) representing the force 
oscillation term is by standard result of type. 


V=A Cos (Wt+ €) 


Where, A = — (1.13) 


ow 
tC 


{(E- «)+ (2}2 


Since the rate at which new notifications actually occur is 
given by yy, we have 





As an alternative to the foregoing approximate discussion Weil) 
we can always investigate special cases by step by step =p (1+vitv) (1.14) 
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Where v1 is the first order solution (1.6). If we take t= 20 
= 68.2 as before thew= 7/20, the putting in (1.13) gives A 
as approximately 8.1p. 


This shows that 10% variation in f', as envisaged by 
Soper, would lead to seasonal variations of about 80% in 
the rate of notifications. 


ALLOWANCE FOR INCUBATION PERIOD: 

Another possible source of error in the use of the simple 
continuous infection model is that is makes no allowance 
for the effect of a fairly well-defined incubation period, 
such as is clearly recognizable is measles (Ludwig, 1974 
and Larger, 1976). 


It we suppose that there is a latent period of length a after 
infection, atthe end of which the infected individual 
becomes infectious, the processes of infection and removal 
then proceeding as before, the deterministic equations will 
be 
dx(t) _ 
dt 
dx(t) _ 
dt 





- Bx(t) y (t-a) + uy (1.15) 





- Bx(t) y (ta) + vy (0) 


The equilibrium vales are exactly the same as for (1.1). If 
we now put D= -, the 
to be. 

(p + ~)u + =~ y=0 


equations for u and v are found 


(1.16) 
-4u+(D+ +-=e-)v=0 
o oc oa 


By shot Heuristic treatment of (1.16), using as operational 
method is as follows. It we assume that the infections 
period is short, both # and yare large and T tends to zero. 


The differential equations for v reduces to 


{(D + =)e% — p} v=0 (1.17) 


The form of the solution depends on the roots of 


(0+ Jer} 


It we write this as 


aD=-log{(p + —)} 


and expand the logarithm in negative powers of D. 


We obtain 


1 1 
aD=-—+ 


oD 202D2 


(1.18) 


The first approximation to (1.18) is given by D = ¥—ao 
and the second 


1 1 
aD=-—+— 
oD 20 


Which is equivalent to the equation 


d?w 1 dw 1 
+——+—w=0 


dt2 2o0dt oa 


(1.19) 


Comparing this shows that the approximate solution here 
is similar to the previous one, expect that we now have a 
instead of t, and the damping coefficient is halved. 


RESULT AND DISCUSSION: 
From the above discussion we see the model has 
represented in special way, in which the infection occurs 


in human body then the resistance f body gradually decays 
same as motion decays in damped oscillation (Kermack, 
and Kendrick, 1927; Hoppenstreadt, 1975). On solving 
the equation of model we get differential equation of 
damped harmonic oscillation, which shows _ the 
phenomenon of occurrence of the infection. Their solution 
gives the idea about the seasonal variation in infection. 
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